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In this paper, the stability of the steady state circular response 
and the bistable response of a flexible rotor — centralised squeeze film 
damper (shortened as SFD) system is analysed, and the system's charac- 
teristics of accelerating or decelerating through the bistable region 
are investigated. It is found that there are two unstable regions for 
the circular response. The larger the unbalance parameter and the 
smaller the bearing parameter and external damping, the easier it is 
to cause an unstable circular response. In addition, the larger mass 
ratio and the smaller stiffness ratio will decrease the threshold rota- 
ting speed of instability. Although in some cases, the system's circu- 
lar response is unstable, the system is still stable and has stable 
nonsynchronous response. 

It is also found that only when the initial rotating speed is 
ahead of the bistable region can the bistable response be produced when 
the system is accelerating through the bistable region. The small 
solution of the bistable response is more stable than the large one. 
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angle acceleration 

bearing parameter (=jjRL 3 /mBu c c 3 ) , where R is the bearing 
radius, L is land length of SFD, and n is the absolute vis- 
cosity of the lubricant 
radial clearance of SFD 
external damping coefficient 
unbalance eccentricity 
retainer spring stiffness 
shaft stiffness 
stiffness ratio (=k a /k s ) 

mass lumped at bearing station and rotor mid-span 
unbalance mass for experiment 
rotating speed in RPM 
unbalance parameter (=e u /c) 
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Subscripts 

(B) = bearing journal 

(D) = disk 

(o) = steady state condition 

(x), (y) = direction subscripts 

INTRODUCTION 

A well designed squeeze film damper (SFD) is a good isolator and 
damper in rotating machinery and can increase the system's stability. 
But poorly designed SFD may destabilize the rotor system, so the sta- 
bility analyses of rotors supported on SFD's have received more atten- 
tion (Ref. l and Ref. 2). The bistable response which is induced by 
the nonlinear fluid film forces will cause a large vibration amplitude 
that is harmful to the rotor system. Since the high nonlinearity of 
the fluid film force makes the mathematical analysis very difficult, 
most stability analyses of the rotor-SFD system are limited to the case 
of steady state circular response of the centralised SFD system (Ref. 3 
and Ref. 4). Whether the system is unstable when its circular response 
is unstable is still an undiscussed problem. As the system is nonlin- 
ear, it may have a stable response other than the circular response. 
Furthermore, it is not clear what the properties are when the rotor 
system accelerates through the bistable region, and the influences of 
the system parameters on the system's stability have not been thor- 
oughly analysed. As these problems are very important for designing 
good SFD, they are worth investigating. 

EIGEN-MATRIX FOR STABILITY ANALYSIS 

For the flexible rotor supported on the centralised SFD shown in 
Fig. l, the nondimensional equations of motion can be derived by taking 
the same assumptions adopted in Ref. 5: 


= Cartesian and polar coordinates 
= x/c, y/c 
= X+iY 

= conjugation of Z 
= mass ratio (=mB/mp) 

= eccentricity ratio of journal or disk (=eB/c or ep/c) , where 
eg and ep are eccentricities 
= damping ratio (=c d /2m I ^ c ) 

- a/»o 

= attitude angles 
= rotor rotating speed 
= pin-pin critical speed (=(k s /mp) 1 / 2 
= speed ratio (=u/« c ) 

= ut for steady state response and w c t for transient response 
= d/dt 
= d/d-c 
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where F=f x +ify is the dimensionless fluid film force. 

In the assumptions of the short bearing and ir-film, there are 
(ref. 6) 
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For the steady state circular response, there are 


i <P 


Zbo =€ o e °i e 0 =const, 0 o =t-<j6 lr 0i=const and Z B0 =iZ B0 


Linearizing the nonlinear fluid film force about the steady state cir- 
cular response Z B0 with the Taylor series, we get (Ref. 6) 

2 i(p — • 2 id) — f 

F=F 0 +A 1 Z 2 +(B 1 +iC 1 )e °Z 2 +A 2 Z 2 + (B 2 +iC 2 ) e °Z 2 (4) 


where, Z2 =Zb-ZbO‘ The linearizing coefficients are 
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As e 0 =const, <j>o=x— 4»i , <|>i=const, it can be seen that these coefficients 
are periodic functions of time, the period being ir. 

In Eq. (1) and Eq. (2), note Z^Z^-Z^, Z 2 =Z B _Z BO' we can 9 et the 


perturbed equations of motion about the steady state response 
Ho DO 

Z BO =e o e and z do = e DCp (approximated to the first order term): 


v 2 0 , 1 

i Zl+ n~ Zl+ n 2 ( z i _z 2 ) =0 
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Thus, the perturbed equations of motion of the system about its cir- 
cular response are in the form of the Hill equation with periodic 
coefficients. 

Now Ui=pi+ip 2 , U 2 =P 3 +ip 4 , and make the transformation of 


Zi=ui ei <*>°, Z 2 =U 2 e i ^ > o (7) 


that is, analysing the problem in the coordinate rotating synchronously 
with the rotor shaft, we have 


10 » > 10 /• k > 1 0 

z l =u l e °, Z 1 =(U 1 +iU 1 )e °, Z ± = (Ui+2iU 1 -U 1 ) e ° (8) 

There are similar results for Z 2 . Substituting Eq. (7) and Eq. (8) 
into Eq. (6) and separating the real and imaginary parts, we get 


P»+DP'+SP=0 (9) 


where P=(Pi,P 2 'P 3 /P 4 ) T ' and the damping matrix D and stiffness matrix S 
are 
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Substituting Eq. (5) into Eq. (9) and noting R=( r n , r 0 , . . . , r Q ) T 

m X .Z O 

= ( P^ ' P^ * P 2 ' Pj> ' P 3 ' P 3 ' P 4 ’ P 4 ) / Eq- (9) can be transferred to the first 
order form: 
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These are the same eight fluid film force coefficients commonly used 
for SFD in the polar coordinate (Ref. 4); then the coefficients matrix 
H is 
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The stability of the system's steady state circular response can 
be determined by the eigenvalue of matrix H. The steady state respon- 
ses Zjjo and Zdo of Eq. (1) and Eq. (2) are asymptotically stable when 
all the eigenvalues of H have negative real parts, and they are 
unstable if only one eigenvalue of H has a positive real part. 

NUMERICAL ANALYSIS OF STABILITY 

The eigenvalues of matrix H are solved with the QR method, and the 
stability of the system's steady state circular response is analysed. 
From numerical results we find that, for both fl<l and Q>1, there are 
two unstable regions. The unstable region generally begins before the 
peak point of response and finishes after the peak point (see Fig. 2). 
In most cases, the larger the unbalance parameter U, and the smaller 
the bearing parameter B and the external damping ratio |3, the easier it 
is to produce the unstable response. Furthermore, the larger the mass 
ratio a and the smaller the stiffness ratio K, the lower the threshold 
rotating speed of instability - perhaps because the larger a and 
smaller K may reduce the critical speed of the system. No instability 
is found in the region of Q<1 and e o <0.5. 

The middle solution of the bistable response is always unstable, 
while the small solution is always stable. For the bistable response 
of the isolated bifurcation form, the large solution is generally 
stable (see Fig. 3), but it may be unstable when the bistable region 
is large (see Fig. 4). For the bistable response of the Duffing form 
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when Q< 1 (although in most cases the large solution is stable) it may 
become unstable in some cases (see Fig. 2). When Q>1, the large solu- 
tion is unstable in most cases (see Fig. 2 and Fig. 3). 

In order to find the reason for this instability, Eq. (1) and 
Eq. (2) are directly solved with the Runge-Kutta method in the unstable 
region. It is found that in most cases, the instability is caused by 
the nonsynchronous response. That is to say, within the unstable 
region, although the circular response is unstable, the system is still 
stable and has a stable, nonsynchronous response. For the example of 
Fig. 3, when Q=2.0, the circular response is unstable, but the system 
has a stable 1/3 subharmonic response at the same time (see Fig. 5), 

Besides the nonsynchronous response, the lack of a steady state 
response is another reason for instability. In only a few cases, the 
system is very unstable, and the response is divergent. 

The above results explain the reason why the bistable response is 
seldom found in practical rotor systems, although it is common in theo- 
retical analyses. 

From the numerical results it is also found that the instability 
can be avoided by selecting reasonable system parameters, especially 
for the instability in Q<1, 

It is worth paying special attention to the fact that the nonsyn- 
chronous response also occurs frequently in centralised SFD, both in 
the regions of Q<1 and Q>1 although the nonsynchronous response may be 
more prevalent in the uncentralised SFD. The nonsynchronous response 
should be avoided as it will cause alternating fatigue to the rotor 
system. 

TRANSIENT CHARACTERISTICS WHEN THIS SYSTEM IS ACCELERATING THROUGH 
THE BISTABLE REGION AND STABLE DEGREE OF BISTABLE RESPONSE 

We have analysed the stability of bistable response. But what is 
the bistable jumping process, and what is its relationship with the 
initial rotating speed? Furthermore, the large and small solutions of 
the bistable response are both stable in most cases, but which is more 
stable? So the transient response of the system accelerating through a 
bistable region is analysed in the following. 

Supposing that the accelerating motion of the flexible 

rotor centralised SFD system has a constant acceleration, that is, 

<{>"= const = a, then the dimensionless equations of motion can be 
written as 

I Z B +2/3Z B +Z D -Z B =U0' 2 e i 0-iUre i 0 

Z B +-i{ (1 +K)Z b -Z d ]=2BF (12) 

2a 
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where x is now denoted as x=co c t and Qg is the initial speed ratio. 

By solving Eq. (12) numerically, we find that, for different ini- 
tial rotating speeds at the beginning of acceleration, the response 
routes when passing through the bistable region are different. If we 
note the beginning and finishing speed ratio of the Duffing bistable 
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response as Q s and fi e (see Fig. 6), then the response will take the 
route of LI when Qq < ^S' and a bistable response occurs (see Fig. 7), 

But the response after jumping down will not completely take the route 
of the steady state response, it goes to steady state after a period 
of oscillation. During the period of oscillation, a nonsynchronous 
response appears in some cases. When the initial speed ratio Qq^S' 
only for very suitable initial response conditions can the response 
take the large solution and cause bistable jumping. In most cases, 
the response will take the route of L2; that is, the small solution is 
the practical one and no bistable jumping occurs. If Qg <fi S' but stops 
accelerating when Q s <<j>Q e / then the bistable jumping down will occur in 
most cases, and when the acceleration is continued, the response will 
take the route of L2. Therefore, the response will take the route of 
L3 for the whole process. 

From the above analysis, it can be seen that the small solution of 
the bistable response is more stable than the large one. In decelera- 
ting, the response takes the route of L4 . It is easy to see that the 
jumping phenomenon of deceleration is not obvious and sudden as that of 
acceleration, and the faster the deceleration, the smaller the jumping 
amplitude . 


EXPERIMENTAL ANALYSIS 

In this experiment, a symmetric rotor is used which has a diameter 
of 32 mm and a length of 889 mm, with a 10.325 kg disk assembled in its 
mid-span. The SFD is fabricated with two load lands, a circumferential 
groove for the oil supply, and no seals. The first rigid critical 
speed of the rotor shaft is about 5500 RPM, the absolute viscosity of 
the lubricant at 18 °C is about ji=8.05xl0 -8 kgs/cm. 

The rotor and journal vibrations in the horizontal and vertical 
directions are measured and recorded by a TEAC-R50C type recorder and 
then processed with an SM-2100B signal analyser. This gives us the 
magnitudes of vibration, the oscillation orbits, and the power spec- 
trum charts of the response. 

Figure 8 is the stable nonsynchronous response of 1/2 the sub- 
harmonic. Although theoretical analyses show that, in this case, the 
steady state circular response of the system is unstable, the system is 
still stable. 

Figure 9 shows a stable synchronous (almost circular) response 
and theoretical analyses also indicate that the circular response is 
stable. 

Figure 10 shows the experimental results of a bistable jumping 
down when system is accelerating through the bistable region. It can 
be seen that the middle solution of the bistable response is unstable, 
as shown in the orbit figure. 


CONCLUSIONS 

1. The perturbed equations of motion of a flexible rotor supported 
on SFD about its circular response are of the form of Hill equations 
with periodic coefficients, and there are two unstable regions for the 
system's circular response separately in Q<1 and Q>1. In most cases. 
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the larger the unbalance parameter and the smaller the bearing parame- 
ter and external damping, the easier it is to produce instability. 
Furthermore, larger mass ratio and small stiffness ratio will decrease 
the threshold rotating speed of instability. 

2. Although the steady state circular response is unstable, the 
system is still stable in most cases and has a stable nonsynchronous 
response. In many cases, especially when Q>1, the large solution of 
bistable response is also unstable, and then bistable jumping will not 
occur. This explains why the bistable response is seldom observed in 
practical rotor systems. Thus further work on bistable response is 
not worth pursuing. 

3. Whether the bistable jumping occurs or not is largely related 
to the initial rotating speed where acceleration begins. Only when the 
initial speed is ahead of the bistable region can the bistable jumping 
be caused. The small solution of the bistable response is more stable 
than the large one. 
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1. squeeze film damper 2. squeeze oil film 

3 . outer race constrained not to rotate 

4. ball bearing 5. centralising spring 

Fig.l Schematic diagram of a flexible rotor 
with squeeze film damper 
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OM 

Fig. 2 Stability analysis of system's circular response 
(U=0.2, £=0.0005, a=0 . 1 , B=0.025, K=0.2) 



'OM 

Fig. 3 Stability analysis of system's circular response 
(U=0 . 4 , £=0.005, < 2=0 . 2 , B=0.05, K=0.01) 
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Fig. 4 Stability of isolated bistable response 
(U=0 . 4 , ;3 = 0.0005, (2=0.01, B=0 . 3 , K=0.001) 


Fig. 5 Nonsynchronous response in 
unstable region (U=0.4, /3=0.005, 
a=0 . 2 , B=0 . 05 , K= 0.01, n=2.0) 


Fig. 6 Schematic diagram of the 
influence of initial rotating 
speed on bistable jumping 
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Fig. 10 System's bistable response by experiment 
(L=1.2cm, c/R=0.1%, m u =4.6g, n=5430rpm) 


62 







